TECHNICAL  PAPER 
NATICK/TP-10/001 


BASIC  PHYSICAL  APPLICATIONS  AND  THE 
MATHEMATICAL  DEVELOPMENT  OF  A  GLINT 
VISUAL  THRESHOLD  DOMAIN  MODEL 


by 

William  J.  Chevalier 


November  2009 


Technical  Paper 
January  1998  -  December  1998 


Approved  for  public  release;  distribution  is  unlimited. 


U.S.  Army  Natick  Soldier  Research,  Development  and  Engineering  Center 

Natick,  Massachusetts  01760-5020 


DISCLAIMERS 


The  findings  contained  in  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army 
position  unless  so  designated  by  other  authorized 
documents. 

Citation  of  trade  names  in  this  report  does  not 
constitute  an  official  endorsement  or  approval  of 
the  use  of  such  items. 


DESTRUCTION  NOTTCF. 

For  Classified  Documents: 

Follow  the  procedures  in  DoD  5200.22 -M,  Industrial 
Security  Manual,  Section  11-19  or  DoD  5200. 1-R, 
Information  Security  Program  Regulation,  Chapter  IX. 

For  Unclassified/Limited  Distribution  Documents: 

Destroy  by  any  method  that  prevents  disclosure  of 


contents  or  reconstruction  of  the  document. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

03-11-2009 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  TYPE 

Technical  Paper 


BASIC  PHYSICAL  APPLICATIONS  AND  THE 
MATHEMATICAL  DEVELOPMENT  OF  A  GLINT  VISUAL 
THRESHOLD  DOMAIN  MODEL 


3.  DATES  COVERED  (From  -  To) 
January  1998  -  December  1998 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

William  J.  Chevalier 


5d.  PROJECT  NUMBER 

85091A 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Natick  Soldier  Research,  Development  and  Engineering  Center 

ATTN:  RDNS-TSM 

Kansas  Street,  Natick,  MA  01760-5020 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


NATICK/TP- 10/001 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


14.  ABSTRACT 

This  paper  discusses  a  theory  that  deals  with  discerning  the  basic  physical  concepts  involved  in  the  attenuation  of  the 
glint  spectral  energy  signature  transmitted  along  two  optical  paths  defined  as  a)  sun  to  ground  reflecting  surface  to 
observer,  and  b)  sun  to  background  reflecting  surface  to  observer.  This  approach  will  provide  a  basis  for  developing  a 
mathematical  model  that  calculates  a  locus  of  points  that  define  a  glint  visual  threshold  domain,  based  on  a  range  of 
specular  surface  reflectivity  coefficients  associated  with  a  corresponding  range  of  solar  incident  angles.  Modifications 
to  the  values  of  the  surface  reflectivity  coefficients  due  to  changes  in  the  optical  properties  of  the  specular  surface  will 
result  in  a  corresponding  change  in  the  glint  threshold  domain  definition.  Thus,  this  model  can  serve  as  a  decision-aid 
tool  for  designing  specular  surfaces. 


15.  SUBJECT  TERMS 

EYE  SUNLIGHT  REFLECTANCE  OPTICAL  MATERIALS 

GLINT  DETECTION  SURVIVABILITY  OPTICAL  PROPERTIES 

MODEL  ALGORITHMS  SKY  BRIGHTNESS  SOLAR 

MODELING  AND  SIMULATION  SPECTRAL  ENERGY  DISTRIBUTION 


SPECULAR  REFLECTION 
SIGNATURE  REDUCTION 
SOLAR  RADIATION 


16.  SECURITY  CLASSIFICATION  OF: _ I  17.  LIMITATION  OF  18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  ABSTRACT  OF  PAGES  Thomas  Gilrov 

U  U  U  SAR  16  19b.  TELEPHONE  NUMBER  (include  area  code) 

508-233-5855 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


UNCLASSIFIED 


This  page  intentionally  left  blank 


UNCLASSIFIED 


Acknowledgments 


The  author  would  like  to  extend  his  appreciation  to  Mr.  Dean-Michael  Sutherland, 
Operations  Research  Analyst,  Science  &  Technology  Directorate,  for  his  timely  and 
constructive  editorial  comments  in  making  this  a  better  paper. 

Special  thanks  are  extended  to  Mr.  David  Tucker,  Operations  Research  Analyst, 
Science  &  Technology  Directorate,  and  Mr.  Bany  De  Cristofano,  Physicist,  Survivability 
Directorate  ,  for  their  constructive  and  timely  technical  comments. 

Special  thanks  are  extended  to  Ms.  Marcia  Lightbody,  our  Chief  Editor,  for  her  very 
constructive  editing  comments  in  helping  to  finalize  this  concept  paper  report. 


v 

UNCLASSIFIED 


This  page  intentionally  left  blank 


vi 

UNCLASSIFIED 


BASIC  PHYSICAL  APPLICATIONS  AND  THE  MATHEMATICAL 
DEVELOPMENT  OF  A  GLINT  THRESHOLD  VISUAL  DOMAIN  MODEL 


1  APPLICATION  THEORY  FOR  MODEL 

The  following  theory  deals  with  discerning  the  basic  physical  concepts  involved  in  the 
attenuation  of  the  glint  spectral  energy  signature  transmitted  along  two  optical  paths 
defined  as  a)  sun  to  ground  reflecting  surface  to  observer  and  b)  sun  to  background 
reflecting  surface  to  observer.  This  approach  will  provide  a  basis  for  developing  a 
mathematical  model  that  calculates  a  locus  of  points  that  define  a  glint  visual  threshold 
domain,  based  on  a  range  of  specular  surface  reflectivity  coefficients  associated  with  a 
corresponding  range  of  solar  incident  angles  (Figures  1  and  2).  Modifications  to  the 
values  of  the  surface  reflectivity  coefficients  due  to  changes  in  the  optical  properties  of 
the  specular  surface  will  result  with  a  corresponding  change  in  the  glint  threshold  domain 
definition.  Thus,  this  model  can  serve  as  a  decision-aid  tool  for  designing  specular 
surfaces. 

The  attenuation  effects  on  transmitted  solar  energy  along  the  optical  paths  within  the 
0.4  -  0.7  pm  visual  wavelength  region  will  be  discussed  in  terms  of  atmospheric 
extinction,  surface  reflective  geometries  and  detector  response.  Reflectivity  properties  of 
a  specular  surface  will  be  presented  as  an  attenuation  function  that  is  wavelength  and 
solar  incidence  angle  dependent.  The  application  of  spectrally  resolved  and  integrated 
concepts,  using  the  contrast  ratio  of  the  glint-to-background  spectral  intensity  values  at 
various  distances  from  the  reflecting  source,  will  be  considered  since  contrast  ratio 
quantification  is  independent  of  solar  variation  intensity. 

The  eye  is  considered  as  the  visual  response  filter  because  we  are  concerned  about 
human  visual  response  to  photo-spectral  stimuli,  although  any  detector  response  filter  can 
be  considered  as  a  response  filter  in  this  model.  Since  psychological  stimuli  and 
responses  are  involved  when  visually  interpreting,  it  is  necessary  to  look  at  the 
photometric  response  of  the  eye  to  the  radiant  flux  density  transmitted  from  a  reflecting 
source  to  the  observer  as  a  spectrally  integrated  weighting  function.  Narrow  band  pass  or 
single  visual  wavelength  filtering  will  be  addressed  as  a  spectrally  resolved  approach  to 
viewing  glint  visual  signatures. 


Figure  1 .  Sun  to  Specular  Surface  and  Background  Optical  Paths 
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Figure  2.  Pictorial  of  Clint  Visual  Threshold  Domains  Based  on 
Two  Sets  of  Specular  Surface  Reflectivity  Coefficients  R(1 )  and 
R(2),  Wavelength  and  Solar  Incidence  Angle  Dependent 


1.1  PHYSICAL  CONCEPT  APPLICATIONS 


1.1.1  RADIOMETRIC  AND  PHOTOMETRIC  DEFINITIONS 

When  undertaking  the  study  of  radiometry  and  photometry,  we  assume  the  existence 
of  an  instrument  called  a  radiometer.  If  radiant  energy  is  incident  upon  a  radiometric 
response  surface  of  known  area  and  orientation  relative  to  its  direction,  then  radiant 
energy  is  converted  to  electrical  energy  in  the  radiometer  as  a  stimulus  response  reaction. 
In  the  same  context,  the  eye  produces  a  bio-metric  stimulus  response  reaction  to  radiant 
energy  in  the  form  of  a  photoelectric  conversion  to  sight. 

Radiometric  quantities  are  physical  quantities  that  are  expressed  in  energy  and 
geometrical  units.  For  purposes  of  this  paper  the  energy  and  geometric  terms  will  be 
defined  as  rate  of  energy  transmitted  per  unit  area,  which  is  called  radiant  flux  density. 
Thus  the  units  can  be  expressed  as  watts/cm2.  The  level  of  radiant  flux  density  will  be 
equated  to  the  level  of  specular  or  glint  reflection  from  a  minor-like  surface,  such  as  eye 
armor. 

The  retina  of  the  eye  is  a  photoelectric  receptor.  Since  perceptual  response  to  physical 
stimuli  is  involved,  the  eye  retinal  receptor  response  to  the  visual  wavelength  spectrum  of 
0.38  -  0.74  pm  could  be  better  defined  as  psycho-physical.  We  are  dealing  with  a  photo¬ 
metric  response  to  a  radiometric  stimuli  of  solar  spectral  energy.  In  this  context,  light  is  a 
visual  aspect  of  radiant  energy,  of  which  the  human  observer  is  aware  through  the  visual 
sensations  that  arise  from  the  stimulation  of  the  retina  of  the  eye.  Brightness  is  defined  as 
that  attribute  of  visual  sensation  by  which  an  observer  is  aware  of  differences  of  observed 
radiant  energy. 

1.1.2  LUMINOUS  ENERGY 

When  the  human  eye  is  used  as  a  photoreceptor  to  the  visual  spectrum  to  measure  the 
relative  levels  of  brightness,  a  relative  luminosity  curve  represented  by  the  function  'V(X)' 
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light  adapted  human  eye.  Generally,  to  convert  the  visual  spectrum  of  radiant  energy  to 
luminous  energy  'Q'  according  to  the  spectral  energy  function  TJjJ,  we  use  the  relative 
luminosity  function  ’V(a,)'  as  a  weighting  function  in  the  following  equations: 

Q=KMV(^)Ul  (1) 

or 

Q  =  KMJo"V(X)Ux<n.  (2) 

where  K.M  is  a  constant  that  determines  the  units  for  Q.  This  equation  provides  the  bridge 
to  convert  radiometric  to  photometric  units.  The  photopic  spectral  luminosity  V(a)  of  the 
human  eye  as  a  function  of  the  wavelength  of  radiant  energy*  is  depicted  in  Figure  3. 

1. 1.2.1  LIMESAL  CONTRAST  RATIOS 

A  highly  developed  sensitivity  of  the  eye  is  its  ability  to  detect  a  small  difference  in 
luminance.  This  difference  is  called  contrast  sensitivity  or  liminal  contrast.  Contrast  'C' 
for  a  given  set  of  conditions  is  defined  as, 

C  =  (Lb  -  L0)/La  (3) 

where  'Lq'  and  'LB'  are  the  object's  or  reflecting  surface's  and  background  luminances 
respectively.  The  luminance  to  which  the  eye  is  adapted  is  'LA'  When  the  reflecting 
object  and  background  fill  the  field  of  view  of  the  eye,  'LA'  is  determined  by  ’L0'  and 
'Lb'.  When  the  reflecting  object  and  background  illuminated  areas  are  approximately 
equal  in  size,  then,  LA  =  1/2(LB  +  L0).  (4) 

If  the  size  of  the  reflecting  or  luminous  object  is  small  compared  to  that  of  the 
background,  which  is  true  regarding  the  eye  armor  reflecting  surface  to  background  size 
ratios,  then  the  eye  adaptation  response  approximates  the  background  illumination 

La=  Lb.  (5) 


*  Erickson,  Ronald,  China  Lake  Report:  "Visual  Detection  of  Targets”,  China  Lake,  1965 
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Wavelength  of  radiant  energy,  um  vs  Luminosity  efficiency 


,  i  i  i  i . i 
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Wavelength  of  radiant  energy,  um 

Figure  3.  Spectral  Response  of  Eye 


1.1.3  EMISSIVITY,  REFLECTANCE,  ABSORPTANCE,  TRANSMITTANCE 

Emission,  reflection,  absorption,  and  transmission  of  radiant  (solar-spectral)  energy 
are  often  treated  as  surface  phenomena.  In  reality,  these  phenomena  exist  when  radiant 
energy  interacts  with  a  surface  of  any  physical  body.  A  true  blackbody,  by  definition, 
absorbs  all  incident  radiant  energy.  However,  a  non-blackbody  is  defined  as  having  an 
absorbtance  ratio  'a(L)'  of  the  absorbed  to  incident  energy  with  a  value  less  than  unity. 
This  means  that  the  balance  of  the  energy  is  reflected  and  transmitted. 

Let  us  assume  that  a  non-blackbody  reflecting  object  having  an  absorptance  ’a(L)'  less 
than  unity,  is  placed  within  an  ideal  blackbody  cavity.  According  to  the  principles  of 
thermodynamics,  the  object  will  reach  the  temperature  T  of  the  cavity  and  remain  at  this 
temperature.  At  this  equilibrium  condition,  the  spectral  irradiance  at  the  object  surface  is 
equal  to  the  spectral  emitted  radiance  'Mj(T)'  derived  from  Planck's  Radiation  Law  as  the 
energy  flow  rate  per  wavelength  per  unit  area  of  a  blackbody  surface.  Yet  the  power 
absorbed  by  the  object  will  equal  'a(X,)M^(T)'.  The  remaining  power  must  be  apportioned 
to  transmission  and  reflection.  If  the  object  is  assumed  opaque,  the  amount  reflected  can 
be  quantified  using  the  expression  [1  -  a(L)]Mx(T).  Because  an  object  is  generally  in 
equilibrium  with  its  surroundings,  it  must  emit  as  much  energy  as  it  absorbs  'a(A.)Mx(T)' 
to  satisfy  KirchhofFs  law.  Emissivity  is  defined  as  the  ratio  of  the  energy  emitted 
by  the  surface  compared  to  the  energy  emitted  by  an  equal  area  of  a  blackbody  surface  at 
the  same  temperature. 

At  the  surface  of  an  object  where  radiant  energy  at  wavelength  T  is  incident  upon  the 
surface,  a  fraction  'a(L)'  is  absorbed,  a  fraction  ’p(a)'  is  reflected,  and  a  fraction  't(a)'  is 
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transmitted.  For  most  materials,  the  radiant  absorbing  value  ’a(A.)'  is  almost  constant  with 
change  in  incidence  angle.  But  the  wavelength  dependent  reflecting  and  transmitting 
radiant  energy  rates  per  unit  area  will  change  with  incidence  angle,  based  on  the  nature  of 
the  optical  characteristics  of  the  materials.  Because  energy  must  be  conserved  we  have: 

a(A,)  +  p(X)  +  x(X)  =  l.  (6) 

The  reflection  may  occur  at  the  surface  by  either  specular  (glint)  or  diffuse  reflection,  or 
it  may  return  from  within  the  material  by  scattering  if  the  material  is  a  translucent,  non- 
homogeneous  medium. 

In  summary,  the  incident,  reflected,  and  emitted  energy  must  be  considered  external  to 
the  object  surface.  Internal  to  the  object  surface,  there  is  absorbed,  transmitted,  and 
scattered  energy. 

1.1.4  ATMOSPHERIC  EXTINCTION  EFFECTS 

There  are  three  atmospheric  processes  responsible  for  the  attenuation  of  transmitted 
optical  images  and  electro-optical  energy  such  as  solar  glint.  They  are  1)  aerosol 
extinction,  2)  molecular  absorption,  and  3)  turbulent  distortion  (scintillation  and  beam 
wander).  Light  propagating  through  the  atmosphere  is  not  only  scattered  and  absorbed  by 
aerosols  and  molecules,  but  the  wave  fronts  are  deflected  and  distorted  by  turbulence. 

The  extinction  has  several  components:  molecular  scattering  and  absorption 
[j(?,)=ji(A.)  +j2(X)]  and  aerosol  scattering  and  absorption  [l(X)=lj(A)  +12(a.)].  The 
extinction  parameter  values  correlate  to  the  loss  of  light  energy  as  it  is  scattered  out  of  the 
beam  or  absorbed  by  the  molecule  and  particulate  constituents  of  the  atmosphere  during 
transmission.  A  combination  of  Bouguefs  and  Beer’s  laws  can  theoretically  quantify  the 
level  of  absorption  of  transmitted  radiant  energy  '<t>(A)’  over  a  distance  rR0'  based  on  the 
nature  of  the  absorption  and  scattering  medium  [k(X)=j(X)+l(A.)]  and  wavelength  of  the 
original  transmitted  energy  '®0(X)'  using  the  following  expression, 

<!>(?,)  =  0>0(X)  exp  [-k(X)Ro].  (7) 

The  distortion  and  tilt  of  image  wave  fronts  by  atmospheric  turbulence  is  represented 
by  *C2n'.  We  can  write  'C2N'  as  a  function  of  temperature  ’C2T'  and  water  vapor  '(C2Q)' 
turbulence  parameters  as  follows: 

Cy  =  (79x1  0*6P/T2)2(C2t  +  0.113  CXQ  +  3.2xlO-3C2Q)  (8) 

where  T’  is  the  pressure  in  millibars,  T  the  absolute  temperature,  and  ’CTQ'  the 
temperature  humidity  co-spectral  structure  function  parameter.  The  refractive  index 
parameter  'Cy  can  be  derived  in  three  ways:  1)  optical  measurement,  2)  measurement 
of  'C2t,  'Ctq'  and  ’C2q',  and  3)  calculation  of  'C2T',  'Cy  and  ’C2Q’  from  bulk 
meteorological  data  made  up  of  water  temperature,  air  temperature,  humidity  and  wind 
speed. 
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The  total  extinction  (a  +  (3)  can  be  measured  optically  by  determining  the  reduction  in 
beam  intensity  over  some  suitable  optical  path.  The  separate  components  can  be 
calculated  from  meteorological  data.  The  molecular  extinction  values  can  be  extracted 
from  a  LOWTRAN  model  and  database  developed  by  the  Air  Force  Geophysics 
Laboratory  (Selby  et  al.,  1978).  The  aerosol  extinction  can  be  calculated  from  the  aerosol 
spectral  density  "N(r  )',  as  follows; 

a  =  /0*  2m2  E(n,X)  N(r )  dr  (9) 

where  r  is  the  particle  radius,  E(n,A),  the  total  scattering  efficiency  at  wavelength  'X'  and 
refractive  index  ’n\ 

1.2  GLINT  CONTRAST  RATIO  MODEL  DEVELOPMENT 

1.2.1  BACKGROUND 

The  application  of  an  imaging  detector,  i.e.,  eye.  Charged  Coupled  Device  (CCD),  etc., 
as  a  device  for  taking  glint  and  background  reflection  measurements  at  various  distances 
from  a  specular  reflecting  source  provides  a  good  field  trial  approach  towards  the 
verification  of  existing  contrast  ratio  models  that  use  spectrally  resolved  and  integrated 
applications.  The  delta  between  experimental  and  model  projected  contrast  ratios  over  a 
distance  from  a  reflecting  source  could  be  significantly  influenced  if  the  model  neglects 
to  consider  the  factor  influences  of  a)  ground  level  atmospheric  extinction  on  solar  energy 
propagation  and  b)  the  effects  of  the  optical  surface's  reflection  coefficients  on  the  level 
of  solar  spectral  energy  reflecting  from  a  surface  as  a  function  of  incidence/reflection 
angle  and  wavelength. 

The  model  assumes  the  glint  reflection  viewed  by  an  observer  is  a  point  source  and  the 
background  an  extended  object.  The  point  source  application  assumes  that  the  spectral 
intensity  decreases  as  an  inverse  square  of  distance.  Glint  to  background  spectrally 
resolved  contrast  ratios  are  independent  of  changes  in  solar  brightness  since  the  glint  and 
background  components  change  simultaneously. 

The  United  Kingdom's  Defence  Evaluation  Research  Agency  developed  a  visual  glint 
contrast  ratio  model  to  calculate  glint  contrast  ratios  over  distances  from  ground 
reflecting  sources  and  uses  these  results  to  feed  into  their  detection  time  model. 

1.2.2  SPECTRALLY  RESOLVED  RELATIVE  INTENSITY 

To  derive  the  representative  expression  for  spectrally  resolved  relative  intensity  we 
need  to  develop  the  relationships  that  describe  the  glint  intensity  viewed  from  a  reflecting 
surface  and  from  a  background. 


*  RI  Young,  RC  Hollins,  T  Holloway  "Simple  Model  for  Predicting  Glint  to  Background 
Contrast  Ratios,  Optical  Glint  Studies  US  and  U.K.  Joint  Venture",  DERA,  U.K.,  July 
1997 
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1 .2.2.1  REFLECTING  SURFACE  TO  OBSERVER 

Assume  the  goggles'  curved  surface  have  radius  'r'  such  that  its  focal  length  is  'r/2'. 

The  goggles  will  reflect  sunlight  into  a  diverging  cone  of  half  angle  (D/2)/(r/2)  =  (D/r), 
where  T)'  equals  cone  diameter  at  the  reflecting  surface.  This  means  that  the  reflected 
cone  at  range  ’R^1  has  a  radius  of:  (R^fTD/r). 

Let  the  glint  reflection  intensity  from  a  goggle  surface  Igog(L)  =  I0(A,,a)Rc(A,,a) 
where:  solar  intensity  function  is  I0(L,a),  and  incident  to  goggle  surface  at  angle  (a), 
and  at  visual  wavelength  (X).  Also,  the  goggle  reflectivity  coefficient  function  is  Rc(A.,a) 
at  visual  wavelength  (X)  and  incident  angle  (a).  Then,  the  reflected  power/unit  area  at 
range  'Rq'  is  defined  as  follows: 

Igog  at  Ro  =  power/unit  area  -  [I0(A.,a)Rc(X,a)7rD2/4  ]/  [tiR^/i2]  =  Igog(7.)r2  /  4R^2  ( 1 0) 

Now  the  sensor  receives  power  over  detector  area  'Ad'  but  lacks  the  angular  resolution 
to  form  a  proper  image  of  glint.  Thus  glint  appears  as  a  feature  of  the  detector  angular 
width  ’9(5’  resolution.  Since  the  perceived  glint  intensity  as  seen  by  the  detector  is 
equivalent  to  the  received  power  divided  by  the  angular  width  squared,  then  the  perceived 
glint  intensity  at  detector  can  be  defined  as  'Igd'.  The  detector  glint  intensity  level  can 
now  be  expressed  as: 

(11) 

which  depends  on  angular  resolution  '@d'  of  the  sensor.  Since  point  source  theory  is 
applied,  the  intensity  relates  inversely  to  the  square  of  the  range  or  distance  'R'  from 
reflecting  source  to  detector. 

1.2.2.2  BACKGROUND  TO  OBSERVER 

The  background  intensity  'Ib(L)'  reflects  with  diffuse  reflectance  'pb  W  times  incident 
solar  intensity  I0(A.)  into  a  hemisphere  over  a  '2tt:'  angle.  From  a  background  of  area  'Ab', 
the  reflected  intensity  at  range  Rq  is 

Ib<X)  at  Ro  =  AhWPbC^  Ro2  =  AbltA)^2  (12) 

The  power  received  by  the  sensor  is  (A^)^)  where  the  perceived  detector  area  is 
defined  by  Ad.  Since  the  background  area  'A^'  is  sufficiently  large  for  a  detector  sensor  to 
form  a  properly  resolved  image  of  angular  size  (A0'5 )/  R^,  then  the  detector  perceived 
background  intensity  Tbd'  can  be  expressed  as  perceived  power  divided  by  the  angular 
width  squared. 

IMW  =  PoXl  AdA„  I„(X)  /  2n  R.1  =  (13) 
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Thus  the  perceived  intensity  of  the  background  by  a  detector  is  independent  of  its 
angular  resolution  and  the  range  between  background  and  reflecting  source,  given  the 
assumption  that  the  diffusion  reflection  has  a  ’2n'  angular  domain  due  to  a  flat 
background. 

1.2.2 3  CONTRAST  RATIO  (GLINT/BACKGROUND) 

We  can  now  define  the  intensity  of  glint  relative  to  the  background  as  perceived  by  the 
detector,  in  terms  of  spectrally  resolved  relative  intensity,  by  solving  for  the  contrast  ratio 
’Cr(a)’  of  the  previously  derived  expressions  of  glint  and  background  reflections  'Igd,  Ibd’ 

CrW  =  MWWW]  -  [WXyibWlOn-2)  /  [2©„: V].  (M) 

This  spectrally  resolved  contrast  ratio  expression  is  independent  of  solar  intensity  but 
depends  on  angular  resolution  of  detector  source,  distance  from  reflecting  source,  and 
optical  reflection  characteristics  of  background  and  reflecting  surface.  The  contrast  ratio 
calculated  over  distance  is  based  on  the  inverse  square  law. 

1.2.3  SPECTRALLY  INTEGRATED  RELATIVE  INTENSITY 

Let  us  consider  the  spectrally  integrated  relative  intensity  as  perceived  by  a  detector 
with  a  spectral  response  function  'f(L)\  The  spectrally  integrated  contrast  ratio  'C{  can  be 
expressed  by  combining  the  spectrally  resolved  expression  'CR(X)'  from  the  previous 
section  with  the  integral  expressions  of  a)  solar  spectral  energy  incident  to  background 
and  glint  surface  T^X,)',  b)  background  and  surface  reflectivity  filter  or  weighting 
functions  'Rfa(L)  and  R^a)’  and  c)  spectral  response  function  of  the  detector  'f(L)’. 

Thus  the  resulting  expression  is 

Cj  =  [ - ](irr2)  /  [2ed2Ro2j.  (16) 

JloMPjA)  f(L)dX 

This  is  a  useful  model  for  ocular  (eye)  detection.  However,  the  model  is  sensitive  to  solar 
spectral  shifts  caused  by  changes  in  the  optical  path  of  sunlight  due  to  the  solar  elevation 
angle  coupled  with  the  atmospheric  extinction  effects. 

Further  analyzing  equation  (16),  the  average  integral  of  I0(L)  =  J  KMV(L)U,dA  is 
equivalent  to  equation  (2).  Based  on  the  application  of  empirical  data  statistical  curve 
fitting  techniques,  the  solar  zenith  and  smooth  surface  reflectivity  coefficient  functions 
highly  correlate  with  the  respective  solar  elevation  and  solar  incidence  to  surface  span  of 
angles  using  a  5th  degree  polynomial  fit  done  in  conjunction  with  the  above  average 
integral  equation.  The  background  reflection  function  'pb(A)'  is  dependent  on  the  nature 
of  the  background  while  the  detector  filter  spectral  response  function  'f(/,)t  is  a  weighting 
function  whose  shape  is  dependent  on  the  nature  of  the  detector  response  filter. 
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1.3  GLINT  VISUAL  DOMAIN  MATHEMATICAL  MODEL 

The  use  of  a  mathematical  model  for  accurately  calculating  a  glint  visual  domain  is 
significant  for  designing  low  reflectivity  optical  surfaces.  This  approach  is  possible  by 
mathematically  describing  the  attenuation  of  solar  glint  energy  over  the  defined  optical 
paths.  An  example  is  the  atmospheric  extinction  effects  on  solar  energy  depicted  by  the 
solar  relative  intensity  versus  zenith  angle  graph*  in  Figure  4. 

ELEVATION  ANGLE  vs  SOLAR  RELATIVE  INTENSITY 


Figure  4:  Atmospheric  Extinction  Effects  on  Solar  Energy 

Another  attenuation  effect  results  from  a  reflecting  surface  interacting  with  incident 
solar  spectral  energy.  A  general  relationship  between  solar  incidence/reflection  angle  and 
surface  reflectivity  coefficients  is  depicted  for  a  basic  eye  armor  surface  in  Figure  5. 

SOLAR  INCIDENCE  ANGLE  vs  REFLECTION  COEFFICIENTS 


CVJ 


SOLAR  INCIDENCE  ANGLE 


Figure  5:  Reflectivity  Characteristics  of  Eye  Armor  Surface 

*Kreith,  Frank,  Principles  of  Solar  Engineering.  Washington,  DC,  McGraw 
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The  final  segment  of  the  optical  path  is  the  atmospheric  extinction  effects  on  solar 
spectral  energy  transmitted  to  an  observer,  based  on  humidity  and  particulate  levels. 
Atmospheric  turbulence  effects  on  transmitted  solar  energy  are  not  considered.  Refer  to 
previous  section  1.1.4:  ’ATMOSPHERIC  EXTINCTION  EFFECTS'  for  details. 

The  spectrally  resolved  and  integrated  contrast  ratio  expressions  developed  in  the 
previous  sections  1. 2.2.3  and  1.2.3  will  be  expanded  to  include  the  effects  of  a)  solar 
zenith  angle  function,  b)  surface  and  background  reflectivity  coefficient  functions,  and  c) 
atmospheric  extinction  absorption  and  scattering  effects. 

1.3.1  SPECTRALLY  RESOLVED  APPROACH 

The  spectrally  resolved  contrast  ratio  expression  can  now  be  expressed  as  follows; 

P0(l)Zc(X,P)Rc(J.,OL)]{exp[-k(>.)R0]}(itr2) 

CRW=[Igd(k)/IM(X)] - - - - - (18) 

PoWPbWPSdV] 

where:  CR(X)  =  contrast  ratio  (wavelength  dependent) 

Igd(A)  =  goggles  reflected  solar  intensity  incident  to  detector 
IM(X)  =  background  reflected  solar  intensity  incident  to  detector 
I0(X)  =  zenith  solar  intensity  incident  to  background  and  reflecting  surface 
Zc(a,P)  =  solar  zenith  angle  coefficient  function  (wavelength  and  zenith 
angle([3)  dependent) 

Rc(X,a)  -  normalized  surface  reflectivity  function  (wavelength  and  incident 
angle(a)  dependent) 

Ph(>-)  =  background  diffuse  solar  reflectivity  function 
I0(X)pb(A.)  =  Ib(X)  reflected  background  solar  intensity 
k(X)  —  atmosphere  total  extinction  coefficient 
Ro  =  distance  or  range  observer  from  glint  reflecting  source 
r  =  radius  of  curvature  of  goggles  (gently  sloping) 

©d  =  angular  resolution  or  spread  of  detecting  source 

By  mathematically  solving  for  Hq'  of  equation  (18),  we  get 

Ro  =  (r/ed){[I0(A)Zc(X,a)Rc(Z,a)/Ib(>.)]{exp[-k(l)R0]}7c  /2C(a)}0'5.  (19) 

By  arbitrarily  setting  values  of  the  following  parameters  as  follows:  a)  the  threshold 
contrast  ratio  ’Cr(a)'  to  1.2,  b)  the  extinction  coefficient  ’k(A)'  to  a  corresponding  level 
of  humidity,  c)  the  solar  elevation  angle  function  'Zc(k,a)'  to  a  solar  elevation  angle,  d) 
the  reflectivity  function  ’Rc(A.,a)'  to  cover  a  range  of  solar  incidence  to  surface  angles,  and 
e)  the  detector  angular  resolution,  a  visual  domain  bounded  by  a  locus  of  visual  threshold 
points  at  ground  distance  Rq  from  the  reflecting  source  which  can  be  generated  by  solving 
for  the  range  Rq.  Empirically  fitted  polynomial  expressions  can  be  used  to  evaluate  the 
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solar  zenith  and  surface  reflectivity  coefficients,  based  on  the  corresponding  solar 
elevation  angle  and  range  of  reflecting  surface  solar  incidence  angles  used. 

Changing  the  reflecting  surface's  optical  characteristics  would  result  in  changing  the 
size  and  shape  of  the  generated  visual  threshold  domain. 

On  perfectly  clear  days  the  extinction  coefficient  *k(X)’  approaches  zero  in  value  such 
that  the  exponential  expression  has  a  value  of  one.  If  there  are  varying  levels  of 
humidity,  then  the  range  Ro  would  have  to  be  solved  iteratively. 

The  advantages  of  using  this  approach  are  that  the  equations  more  accurately  describe 
the  basic  physical  phenomena  of  the  optical  path. 

As  with  the  application  of  any  contrast  ratio  concept,  the  field  trial  verifications  would 
be  generally  independent  of  the  energy  intensity  level  fluctuations  but  depend  on  detector 
source’s  angular  resolution.  The  intensity  decreases  as  the  inverse  square  of  range. 
Regarding  the  reflecting  surface  geometry,  the  reflecting  surface  effective  radius  of 
curvature  is  the  square  root  of  the  product  of  the  radii  components. 

1.3.2  SPECTRALLY  INTEGRATIVE  APPROACH 

We  can  use  the  spectrally  integrative  approach  even  though  the  contrast  ratio  is 
sensitive  to  spectral  changes  such  as  diurnal  spectral  shifts.  The  expression  would  be  an 
integral  over  the  0.4  -  0.7pm  eye  photo  visual  region  of  the  solar  spectrum  of  the  previous 
spectrally  resolved  expression. 

It  is  expressed  as; 

fl0(A)Zc(X,(3)Rx;(X,a)f(A)d(X) 

q  =  [- - - ]{exp[-k(X)Ro  ]}(m2)  /  [20^]  (20) 

Jl0(Z)Zc(X,P)pb(l)f(X)d(Z) 

where:  I0(A.)  =  zenith  solar  intensity  (wavelength  dependent) 

Zc(a,P)  =  solar  elevation  angle  function  (wavelength  and  zenith  angle 
dependent) 

Rc(A.,a)  =  surface  reflectivity  function  (wavelength  and  incident  angle 
dependent) 

p(L)  =  background  diffuse  reflectivity  function 

k(X)  =  atmosphere  total  extinction  coefficient 

Ro  -  distance  or  range  observer  from  glint  reflecting  source 

r  =  radius  of  curvature  of  goggles  (gently  sloping) 

9d  =  angular  resolution  or  spread  of  detecting  source 
f(X)  -  detector  spectral  response  function 
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We  solve  the  spectrally  integrated  expression  equation  (20)  for  threshold  range  or 
distance  Rg  from  the  reflecting  source  as  we  did  with  the  spectrally  resolved  equation 
(18). 

Values  are  set  for  the  parameters  like  we  did  when  using  the  spectrally  resolved 
equation  as  follows:  a)  the  threshold  contrast  ratio  ’CR(A.)'  to  12,  b)  the  extinction 
coefficient  ’k(A.)'  to  a  corresponding  level  of  humidity,  c)  the  solar  elevation  angle 
function  ’Zc(X,a)'  to  a  solar  elevation  angle,  d)  the  reflectivity  function  'R^cc)'  to  cover 
a  range  of  solar  incidence  to  surface  angles,  and  e)  the  detector  angular  resolution. 

Regarding  optical  path  evaluation,  we  integrally  average  the  incident  solar  spectral 
energy  function  at  zenith  to  obtain  the  solar  luminosity  incidence  value,  i.e.,  140,000 
candles/cm  .  We  evaluate  the  solar  zenith  angle  and  surface  reflectivity  coefficients 
using  the  corresponding  known  and  'empirically  fitted'  polynomial  equations. 

The  product  of  the  solar  zenith  incidence  energy  ’I0(?-)’ ,  solar  zenith  angle  coefficient 
'Zc(A,p)'  and  the  range  of  surface  reflectivity  coefficients  rRc(A,a)'  yield  a  range  of  glint 
surface  reflection  values  'Igog'  that  are  further  attenuated  by  the  atmospheric  extinction 
effects  on  the  glint  signature  propagated  over  the  distance  to  observer.  The  filter 
weighting  function  'f(X)'  and  the  angular  resolution  '0d'  of  the  detector  determine  the 
detector  response  to  the  incoming  glint  and  background  energy. 

The  background  spectral  energy  incoming  to  the  detecting  source  is  the  product  of  the 
incident  to  background  surface  solar  spectral  energy  10(X)ZC(A,,|3)'  and  the  background 
diffuse  reflectivity  function  'p(X)'  attenuated  by  the  atmospheric  extinction  effects. 

A  visual  domain  bounded  by  a  locus  of  visual  threshold  points  at  ground  distance  Rg 
from  reflecting  source  is  generated. 

1.3.3  MODEL  APPLICATIONS 

This  model  can  be  converted  to  a  software  package  and  used  as  an  optical  surface 
design  tool  by  allowing  the  user  to  evaluate  how  changes  in  the  optical  characteristics  of 
reflecting  surfaces  will  impact  on  the  size  and  shape  of  the  glint  visual  threshold  domains. 

The  model  accurately  interprets  the  attenuation  of  transmitted  solar  energy  along  its 
optical  path,  thus  enhancing  realism  as  a  decision  aid  tool  for  surface  design. 

1.3.4  FUTURE  CONSIDERATIONS 

The  effects  of  atmospheric  turbulence  on  glint  will  be  evaluated  in  the  future  by 
modeling  gaussian  shaped  glint  noise  scintillation  modulated  to  a  continuous  wave  glint 
signature.  There  is  a  potential  need  to  model  and  verify  significant  differences  in  the 
observer's  signature-acquired  detection  times  at  a  reasonable  probability  of  detection  level 
between  gaussian  and  nongaussian  glint  signatures  at  various  observer  distances  from 
reflecting  source.  For  example,  a  significant  difference  surface  threshold  domain  can  be 
generated  from  a  field-generated  data  base,  using  statistical  fitting  techniques,  that  depicts 
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statistically  viable  gaussian  and  nongaussian  detection  time  differences  against  glint 
intensity  and  glint  signal  to  gaussian  spectral  noise  levels. 

A  key  question  that  must  be  answered  is  "under  what  combat  scenerio  circumstances 
does  glint  hazard  become  significant  enough  to  compromise  the  mission"?  If  there  are 
circumstances,  "what  range  of  significantly  different  gaussian  minus  nongaussian  glint 
values  located  within  the  threshold  contour  domain  would  contribute  to  the  glint  hazard" 
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